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c 2 ^ , Abstract. We consider the continuous time, zero-temperature heat-bath dynamics for the 

04 ■ nearest-neighbor Ising model on 1?" with positive magnetic field. For a system of size _L £ N, we 

start with initial condition a such that Ox = —1 if a; € \—L^V^ and fJx = -1-1 and investigate 
■ the scaling limit of the set of — spins when both time and space are rescaled by L. We compare 

\^ ' the obtained result and its proof with the case of zero-magnetic fields, for which a scaling result 

was proved in [15]. In that case, the time-scaling is diffusive and the scaling limit is given by 
anisotropic motion by curvature. 



Oh; 

■ 1. Introduction 

> ■ 

J2 ■ The Ising model is one of the simplest model proposed by statistical mechanics to investigate 

ferromagnetic properties of metals. It is based on the following simplification of reality: we 
assume that a ferromagnetic material is composed of macroscopic magnet that live on a lattice 
and can assume only two orientation, up (or -|-) and down (or — ) called spins; the way orientation 
^ . of the micromagnet are determined follows two rules: neighboring magnets like to have the same 

orientation, and all spins like to align with the external magnetic fields if there is one. 

The Ising model is a probabilistic model (i.e. a probability law on the set of possible spin 
orientation) following these rules and is defined using a Boltzman-Gibbs formalism, which is 
detailed in the next section. It describes the equilibrium state of a ferromagnet. 

CN ■ The Stochastic Ising model also called Glauber dynamics or heat-bath dynamics for Ising 

model, is a Markov chain on the set of spin configuration that describes the evolution of a 
magnet out of equilibrium (e.g. after a brutal change of temperature or of external magnetic 
^ . fields). This is the object we study in this paper. Our aim is to understand the time needed 

^ \ and the pattern used for a magnetic material to reach its equilibrium state after a change of 

external conditions. In order to bring a more complete answer to these questions, we consider 
them in a simplified but still non-trivial setup that is the zero-temperature limit. This brings 
heuristic intuition for what should happen in the whole low temperature regime. 
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Z+i :={x + (l/2) I xeZ}. 



1.1. The Glauber Dynamics for the Ising Model. Set 

Z* 

Consider the square 

A = Al := [-L,L]2n(Z*)2. 

We define the set of spin configuration on A/, to be { — 1,1}"^^. A generic spin configuration 
is denoted by o" = {(yx)x^Ki, and Ox G {— 1} is called the spin at site x. Define the external 
boundary of A as 

aA:={yGZ'^\A I 3xG A,x~y}. (1.1) 
The Ising model at inverse temperature /3, with external magnetic fields h and boundary 
condition 77 G {— Ij 1}^"^^ is a measure on the set of spin configurations defined by 
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and the convention is taken that ax ■= ijx when x ^ A. Note that the first term in the exponential 
makes the spins of neighboring sites more hkely to agree whereas the second term underlines 
that configuration with spins aligned with the magnetic fields are favored. 

The heat-bath dynamic for the Ising model (at inverse temperature /?, with external magnetic 
fields h and boundary condition rj) is a Markov chain on the set of spin configuration { — 1, 1}"^. 
We denote the trajectory of the Markov chain by a{t) = {(Txi't))xGA- One starts from a given 
configuration ao and the rules for the evolution are the following: 

(i) Sites X £ A are equipped with independent rate one Poisson process: (T^)n>o where 
Tq = 0, i.e. the increments {Tn^x — '''n-i,x)n>o,x&<^ exponential variables. 

(ii) The spin at site x may change its value only when the clock at x rings, i.e. at time r^, 
71. > 1, and its value stays constant on the intervals \Tn^'^n+i)^ ^ 

> 0. 

(iii) When a clock rings, the spin at x is updated according to the following law independently 
of the past history of the process: 

— cra:(io) = + with probability 



<7x(io) 



exp (pT.y^x^y{'^n,x) + 
2 COsh(/3 Yjyr^x ^yi'^n,x) + h) 
with probability 

exp {-p Ey~x c^j/lVx) + h 



(1.4) 



2 COsh(/3 Y^yr^x ^^S/(Tn,x) + h) 



(1.5) 



When X has a neighbor y in dA, the value ay appearing in the equations (1.4)-(1.5) is fixed 
(by convention) to be equal to r]y for all time, where rj is the boundary condition. The update of 
the spin corresponds to sampling a spin configuration according to the measure fi^^'^{- \ ay = 



As a consequence to this last remark, fj,^^'^ is the unique invariant measure for the dynamics, 
so that the law of a{t) converges to the equilibrium measure 

The main questions in the study of dynamics is how much time is needed to reach equilibrium 
and what is the pattern used by the chain to reach it when we consider the dynamics on a very 
large domain A . The answer to this question should of course depend on the temperature /3, 
the magnetic fields h, the boundary condition rj the and the initial condition aQ. 



In what follows, we denote often by + resp. 
where all spins are + resp. — . 



the spin configuration or boundary condition 
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1.2. Conjecture and known results in the case /3 S (0,oo), h = 0. Let us review shortly 
what is known and conjectured for these kind of dynamics when /3 G (0, oo) and h = 0. The 
property of the dynamics depends crucially on the equilibrium property of the system. Recall 
that the two dimensional Ising model undertakes a phase transition at /3c = log(l + V^)/2 (see 
the seminal work of Onsager [22]) which has the following form: 

• When P < f3c, correlation between spin at different sites decay exponentially, and for 
this reason what happens in the center of becomes independent of the boundary 
condition when L tends to infinity. This is called the high temperature phase. 

• When /3 > f3c, on the contrary, long-range correlation are present between spins, and 
boundary condition plays a crucial role. In particular the measure fJ-^^'^ and /i^''*'" 
corresponding to + and — boundary condition are very different. This is called the low 
temperature phase. 

In the high temperature phase, the rapid decay of correlation between distant sites makes the 
evolution of the system in two distant zones of the box Al almost independent. Schematically, 
the box of Al can be separated in 0{L^) zones of finite size that come to equilibrium indepen- 
dently. This requires a time of order log L. This prediction has been made rigorous by Lubetzky 
and Sly, who proved that in that case, the mixing time (i.e. the time to reach equilibrium) for 
the dynamic in Al is equal to Aoo logL(l + o(l)) where Aqo is the relaxation time for the infinite 
volume dynamics (see [18]). For a formal definition of mixing time and relaxation time and an 
introduction to the modern theory of Markov chain, we refer to [16]. 

In the low temperature phase the behavior of the dynamics depends on the boundary condition 
and for the sake of simplicity we restrict to the case of + boundary condition. In that case, the 
equilibrium state is biased towards plus, and even spins in the center have a larger probability 
to be + than minus. In a sense, one can say that, at equilibrium, the center of the box "knows" 
what the boundary condition is. Thus, if one starts e.g. from full — initial condition (crx(O) = — 1 
for all X £ A), information must travel from the boundary to the center of the box in order to 
reach equilibrium. For this reason the mixing time is much longer than in the high temperature 
phase. 

In [17], Lifshitz described a conjectural pattern used by the system with + boundary condition 
to reach equilibrium that can be described as follows: starting from — initial condition the system 
should rapidly reach a state of local equilibrium that looks like the equilibrium measure with — 
boundary condition, (we call this the — phase); then on the time-scale L^, something looking like 
the true equilibrium measure with + boundary condition, the + phase, should start to appear in 
the neighborhood of the cubes boundary. The interface between the + and the — phase should 
move on the diffusive time-scale L^, having a drift in time proportional to its local curvature. 
As a consequence the system should reach equilibrium when the bubble formed by the — phase 
disappears macroscopically, i.e. in a time 0{L?'). 

For finite (3 > (3^ this conjecture is far from being on rigorous mathematical ground, but 
Lifshitz ideas have been used to get bounds on the mixing time. The best to date being by 
Lubetzky et al. [21] saying that the system reaches equilibrium in a time L}°^^, still far from 
the conjecture L^. This gap between the Lifshitz conjecture and the rigorous mathematical result 
has been one of the incentives to study the simpler zero-temperature version of the model. 

Dynamics with different boundary condition or non-zero magnetic fields at low temperature 
also exhibits interesting behavior like low-temperature induced metastability (see e.g. [26]) that 
we choose not to expose here. Note also that results are available at the critical temperature /3c 
[19] , where the equilibrium state of the system is somehow harder to describe. 
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2. Zero-temperature dynamics 

Given h > 0, we look at the hmiting dynamics when /3 tends to infinity. We call this the 
zero-temperature limit (recall that (5 is the inverse of the temperature). In what follows we will 
consider only -|- boundary condition: rjx = +1, Vx G dA. 

With this setup, on a finite box with -|- boundary condition, the limit of the Ising measure 
lim/3_j.oo /i^'^'"'' is just the Dirac measure on the full -|- configuration: ax = +1, Vx E Al. The 
limiting dynamics is non-trivial and can be described as follows: the value the spin at x is 
updated with rate one as before; when a spin is updated it takes the value of the majority of 
its neighbors if it is well defined and take value ±1 with probability e^'*/(2 cosh(/i)) if it has 
the same number of + and — neighbors. The process can also be defined when h = oo, it 
corresponds to the case where both (3 and h tend to infinity with h <^ (3. 

The question concerning the pattern used to reach the equilibrium then takes the following 
form: starting from a finite domain of I? filled with — spins what is the time needed to reach 
the whole + configuration and what is the pattern used to reach it. 

More precisely: We consider a compact, simply connected subset V C [—1, 1]^ whose boundary 
is a closed smooth curve. Given L € N we consider the Markov chain described above with initial 
condition 

..(0) = (-> (2.1, 

I +1 otherwise. 

In order to see a set of "— " spins as a subset of M^, each vertex x E (Z*)^ may be identified 
with the closed square of side-length one centered at x, 

Cx:=x + [-l/2,l/2]\ (2.2) 

One defines 

AL{t):= U Cx, (2.3) 

{x: a^{t) = -l} 

which is the "— droplet" at time t for the dynamics. Note that the boundary of Aiit) is a union 
of edges of (and this is the reason why we defined the Ising model on (Z*)^ rather than on 
Z2). 

We investigate the scaling limit of ALit) and the time needed for it to vanish. The nature of 
the scaling limit depends really much on the value of the external magnetic fields h: 

• If /i > 0, the interface between + and — is always pushed towards the minus region, at 
linear speed, and thus macroscopic motion is visible on time-scale L. 

• If /i = 0, there is no sign that is favored, and the interface will be pushed in the direction 
of its curvature to reduce its length. This is visible only on the diffusive time scale. 

The main aims of this paper are: to review in detail the recent proof (with co-authors [15]) 
that the scaling limit of Aiit) is given by an anisotropic curve shortening fiow, and to give a 
description of the scaling limit in the case h > (for the sake of simplicity we will limit our 
proof to h = oo). 

2.1. The case of zero magnetic fields. Let us focus first on the case of zero-magnetic fields, 
and make more precise the Lifshitz conjecture [17] concerning low temperature dynamics in this 
case. 

On heuristic grounds, Lifshitz predicted that the boundary of ^^(t) should follow an anisotropic 
curve shortening motion: after rescaling space by L and time by and letting L tend to in- 
finity, the motion of the interface between ALit) and its complement {i.e. between -|- and — 
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spins) should be deterministic and the local drift of the interface should be proportional to the 
curvature. An anisotropic factor should appear in front of the curvature to reflect anisotropy 
of the lattice. Spohn [27] made this conjecture more precise and brought some elements for its 
proof: let "y{t,L) denote the boundary of the (random) set {1/ L)AL{L'^t). Then, for L — )■ oo, 
the flow of curve {'y{t, L))t>o should converge to a deterministic flow {'y{t))t>o and the motion 
of the limiting curve should be such that the normal velocity at a point x G 7(t) is given by the 
curvature at x multiplied an anisotropic factor a{9x), where 9x is the angle made by the outward 
directed normal to 7(t) at x together the horizontal axis (see Figure 1). The velocity points in 
the direction of concavity. The function a(-) has the explicit expression 

"^^^ 2(|cos(e)| + |sin(0)|)2- ^^-^^ 

The area enclosed by the curve decreases with a constant speed: — f^^ a{0)d0. In particular, 
the curve j{t) shrinks to a point in a finite time 

_ AreaiV) _ Area{T>) 

Note that the function a(-) is symmetric around and is periodic with period 7r/2. This is 
inherited from the discrete symmetries of the lattice (Z*)^. 




Figure 1. A graphic description of the support function h. Given 6, consider 
the point x{9) of 7 that maximizes x ■ v{6) (it is unique if the curve is strictly 
convex). Then h{0) = x{6) ■ v{9), and k{9) is the norm of the curvature vector of 
7 (green vector) at x{9). If the tangent to 7 at x exists it is normal to v{9) and 
\h{0)\ is the distance between the tangent and the origin. 



In a recent work [15] , with co-authors, this conjecture was brought on rigorous ground in the 
case when the initial condition D is convex. Let us first give a rigorous definition of motion by 
curvature and a result concerning its existence: given a convex smooth closed curve 7 = dT> 
in M^, we parameterize it following a standard convention of convex geometry (cf. e.g. [8] and 
Figure 1). For 9 G [0, 27r] let v{9) be the unit vector forming an anticlockwise angle 9 with the 
horizontal axis and let 

h{9) = sup{x ■ V (9), xej}, (2.5) 
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where • denotes the usual scalar product in M?. The function 9 i— t- h{6) (called "the support 
function") uniquely determines 7: 

V= Pi {x : X -vie) < h{e)}. (2.6) 

0<e<27r 

With this parameterization, the anisotropic curve shortening evolution reads 

(dth{e,t) = -aie)k{e,t), 
\h{e,o) = h{9), 

where, for a convex curve 7, k{6) > is the curvature at the point x{6) £ 7 where the outward 
normal vector makes an angle 9 with the horizontal axis and the time derivative is taken at 
constant 9 (see [8, Lemma 2.1] for a proof that (2.7) is equivalent to the standard definition of 
anisotropic motion by curvature). 

Existence of a solution to is not straight-forward. It was proved under the assumption that 
a is in [8]. The function a given by equation (2.4) is only Lifshitz due to singularity at 
9 = 111/2, i = 1, . . . ,4 but a proof of existence and uniqueness of the motion was given in [15] 
for that particular case. 

Theorem 2.1 ([15, Theorem 2.1]). Let V C [—1, 1]^ he strictly convex and assume that its 
boundary 7 = dD is a curve whose curvature [0, 27r] 9 ^ 1— )• k(9) defines a positive, 2Tr-periodic, 
Lipschitz function. 

Then there exists a unique flow of convex curves {'j{t))t with curvature defined everywhere, 
such that 7(0) = 7 and that the corresponding support function h{9,t) solves (2.7) fort > and 
satisfies the correct initial condition h{9,0) = h{9). 

The curve 7(t) shrinks to a point xj G at time tf = Area{T>)/2. 

For t < tf, 7(t) is a smooth curve in the following sense: its curvature function k{-,t) is 
Lipschitz and bounded away from and infinity on any compact subset of [0,tf). 

Now let Vit) denotes the flow of convex shape whose support function is solution of (2.7), 
with T>{0) = T>. Let B{x, r) denote the closed Euclidean ball of radius r > and center x G M^. 
For a closed set C C , define the inner and outer (5- neighborhood of C to be 



■.= \jB{x,5), 



( \ - (2.8) 

We say that an event, or rather, that a sequence of events (^n)n>o holds with high prob- 
ability of w.h.p. if the probability of tends to one. Then P(t) is the scaling limit Ah^t) in 
the following sense. 

Theorem 2.2 ([15, Theorem 2.2]). Consider the dynamics starting with initial condition given 
by (2.1), where D is a convex shape satisfying the assumption of Theorem 2.1. For any (5 > 
one has w.h.p. 

V^-^\t) C yAiiLh) cV^^\t) for every 0<t<tf + d (2.9) 
AiiL^t) = for every t>tf + 5. (2.10) 
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In particular, one has the following convergence in probability: 

lim ^ , , = -. (2.11) 

We can mention previous related result in the literature: in [4], the authors consider a simpli- 
fied dynamics that do not allow the interface to break in several components, for this dynamics, 
they present a result similar to (2.11) without any statement concerning the limiting shape. 

In [3] the drift for the interface at the initial time is studied, and the author prove that it is 
proportional to the curvature multiplied by an anisotropy function that is different from a{6) of 
equation (2.4). This difference is explained by the fact that the initial condition considered by 
the authors (2.1) is very far from being a local equilibrium for the interface dynamics. 

Remark 2.3. Our result concerns only the case where P is a convex shape and there are serious 
reason for it. The first one is that there is no proof of the existence of anisotropic curvature 
motion starting from non-convex initial condition (see [9] in the isotropic case) , and its regularity. 
On the probabilistic side, our method seem to be able to handle the non-convex case, but would 
need considerable improvements. 

We can now compare this result with what happens in the case of positive magnetic fields. 

2.2. Zero-temperature dynamics w^ith positive magnetic fields (/i = oo). Consider the 
dynamics where sites are still updated with rate 1, but with the following rule of update: when 
a site is updated, its spin flips to -|- if it has two or more -|- neighbors and to — if has three or 
more — neighbors. 

In that case, the right time-scale to describe the evolution of is not but L and the 
interface is always going towards the — side regardless of the curvature. The main new result 
of this paper is identifying the scaling limit in this case. 

Heuristically, the intensity of the drift of the interface can be deduced from mathematical 
work on Totally Asymmetric Simple Exclusion Process (TASEP), more precisely from results 
giving the scaling limit of the height function (see Section 3.2 for detailed explanations). The 
drift at a point where the interface makes an angle 6 with the horizontal axis is equal to 

,fn^. |sin(2g)|(|cosg| + |sing|) 

^^^^ ■= l + |sin2e| • ^^-^^^ 

Let us give a rigorous definition of this shape evolution in the case of convex initial condition. 
The shape remains convex at all times and one can describe the evolution of the interface in 
terms of the support function (recall (2.5)) as follows: 

(dth{e,t) = -b{6), 

\h{e,o) = h{e). ^ ■ ' 

The problem is that the equation (2.13) is not well-posed. However there is a notion of weak 
solution to (2.13) that has a rather simple description and this is the one that will be of interest 
to us: given an initial condition P, define 

V{t) := Pi {x G ]R2 . ^ . ^(0) < _ (2.14) 

eG[0,27r] 

It could be shown that the support function of {'D{t))t>o is the unique viscosity solution of 
(2.13) (see [6] for an introduction to this concept) but we will not pain ourselves with such 
considerations, and rather consider (2.14) as a definition. 
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Our main result is 

Theorem 2.4. The process {j^AL{Lt))t>o converges to {'D{t))t>o probability, in the topology 
of time-uniform convergence for the Hausdorf metric. This is to say, for any 5 > with 
probability tending to one, 

V^-'\t) C (^^L(Lt))i>o C V^'\t). (2.15) 

For simplicity we choose to expose the proof in the case where the initial condition is the 
square [—1, 1]^. Starting from a general convex domain is not conceptually more complicated 
but would involve notational complication that we wish to avoid. 

Let us describe T>{t) when P = [—1, 1]^: define the function 51 on M x M_|_ by 

:= if 1^1 (2.16) 

I |x| if |x| > t, 

let Vi{t) denote the epigraph of g(-,t) - 2 in the base (0,fi,f2) where fi := h ■= ^"4^' 

i.e. 

Vi{t) :={xfi + yf2 I xeR,y>g{x,t)-2}, (2.17) 
and T>i, i = 2, . . . , 4, denote its image by rotation of an angle {i — l)vr/2. Then one has 

4 

P(t) := f|A(t). (2.18) 
i=l 

It is a convex compact set and it is the solution of (2.14) when h{9) = \/2/(| cos(6' + 7r/4)| + 
I sin(6' + 7r/4)|) is the support function of [—1, 1]^. 

Remark 2.5. The proof of Theorem 2.4 also adapts quite easily to the case of positive magnetic 
fields h G (0,00). In that case the scaling limits remains the same but time has to be rescaled 
by a factor cosh(/i)/sinh(/i). The non-convex case is more delicate as the limiting shape has not 
a nice description and it might split into several connected components. All of this makes the 
analysis more complicated on a technical point of view but we believe that an analogous result 
could be proved in that case with some efforts if the initial condition is sufficiently nice, e.g. 
with smooth boundary. 

2.3. Interpolating between h = and h > 0: the weak magnetic field limit. Remark 
2.5 says that the scaling limit of set of — spins is somehow independent of the intensity of the 
magnetic fields apart from an armless scaling factor. We want to discuss shortly here a way to 
obtain a non trivial intermediate regime (sometimes referred to as a crossover regime) between 
the cases h > and h = 0. We give a brief description of a conjecture concerning that case, 
based on heuristic consideration. 

The intermediate regime should take place when h = Hl ^ aL~^ . When h ^ Theorem 
2.4 should hold whereas Theorem 2.2 should be valid when h <^ . In this intermediate 
regime the scaling should be LP' and the limiting equation for the support function 

{dth{e, t) = -a{0)k{e, t) - ab{e) 

\hi9,0) = h{e) ^ ' ^ 

The existence and uniqueness of solution to the above equation is a challenging issue. 
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Note that this type of scaUng has already be studied in the context particle system in the 
case where the interface between + and — is the graph of a function (see Section 3.2) and is 
often referred to as Weakly Asymmetric Exclusion Process or WASEP. 

It has been proved in [5, 10, 12] that in a certain sense, the height function of the particle 
system converges to the solution of 



which is as can be checked by the reader, is equivalent to (2.19) in this case. These result could 
be used as a starting point to get the scaling limit starting from a compact convex shape filled 
with — . Again, the major difficulty seem to be on the analytical side, as on the probabilistic 
side, most of the tools applied in the proof of Theorem 2.2 could be applied. 

2.4. Higher dimensions. Of course, the Ising model and the dynamics are also defined in 
higher dimensions d > 3, for which one as also a phase transition, and they have also raised a 
lot of interest (in particular the case d = 3 for obvious reason). Let us briefly talk about what 
should be true in that case with a focus on low-temperature with plus boundary condition (for 
high temperature, we can notice that the result of [18] holds in every dimension). 

The Lifshitz conjecture for the low-temperature dynamics should hold in fact in any dimension. 
However rigorous results are even more difficult to obtain for the model with zero-magnetic fields: 
for finite /3 the best known bound for the mixing time at low temperature with + boundary 
condition is super-exponential (it is equal to 0(exp(L'^~^(log L)^)) [28]), which is really far from 
the conjectured 0{L?') 

Concerning the zero-temperature case, results are much more precise: in [2], it has been proved 
that for d = 3, the time necessary for the last — spin to disappear starting from a full cube 
with + boundary condition is of order L^(logL)*^(^), a similar upper bound has been derived 
for arbitrary dimension d in [13] (with no matching lower bound). However, it seems that with 
actual tool, we are very far from being able to prove a shape theorem similar to Theorem 2.2 
when d = 3, or even to get rid of the log. A first step would be to derive on a heuristic level the 
anisotropy function (the equivalent of a[9) from (2.4)), and it seems highly non-trivial. 

The anisotropy function should have a lot of singularity when d > 3. For instance, consider 
the three-dimensional dynamics on (Z*)^ starting from initial condition 



where .6(0, L) denote the two dimensional Euclidean ball. If L > 4, then this initial configuration 
is stable, and the system stays forever in that state: every — spin has a strict majority of agreeing 
neighbors and the same for + spin. However, the mean curvature is well defined and positive at 
every point of the interface of the cylinder. 

For the case of zero-temperature with positive h, it should be rather simple to show that an 
hyper-cube full of minus of diameter L needs a time of order L to disappear. Getting the exact 
asymptotic and a shape Theorem seems more difficult but probably not out of reach. 




(2.20) 



a{x) = — 1 if X G 
(t{x) = +1 else , 




(2.21) 
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2.5. Organization of the paper. In Section 3 we present some of the main tools of the proof. 
They are: general monotonicity property of the dynamics, correspondence with particle systems 
and scaling limits for interface dynamics. 

In Section 4, we sketch in detail the proof of Theorem 2.2 from [15] 

In Section 5, we prove our main result, that is Theorem 2.4, and underline the similarities 
and difference between this proof and the one of the zero-magnetic fields case. 

3. Interface dynamics, correspondence with particle systems and other 

technical tools 

3.1. Graphical construction and monotonicity. We present here a construction of the 
dynamics (called sometimes the graphical construction) that yields nice monotonicity properties. 
We consider a family of independent Poisson clock processes {t^)x£Z2. that is, to each site 
X € (Z*)^ one associates an independent random sequence of times {Tn)n>o, that are such that 
Tq = and (t^_|_i — T^)n>o are HD exponential variables with mean one. One also defines 
random variables (f/n,a;)n>Oi 2; G (Z*)^ that are IID Bernoulli variables, that assume value ±1 
with probability e^^ / {2 cosh.{h)) . 

Then given an initial configuration ^ G { — 1, 1}^^*^^ one constructs the dynamics cr^(t) starting 
from (7^(0) = as follows 

• i(^x{t))t>o is constant on the intervals of the type [t^,t^_|_i). 

• ax{Tn) is chosen to be equal to ±1 if a strict majority of the neighbors of x satisfies 
'^yi'^n) — and Un,x otherwise (this definition makes sense as, almost surely, two 
neighbors will not update at the same time). 

This construction gives a simple way to define simultaneously the dynamics for all initial 
conditions and boundary condition, using the same variable U and r (we denote by P the 
associated probability). Moreover the coupling of dynamics thus obtained preserves the natural 
order on { — 1, +1}^^*)^, given by 

C>C'^^x>Cx for every x G (3.1) 

(this order is just the opposite of the inclusion order for the set of "— " spins, which is therefore 
also preserved). Indeed, if ^ > ^' and resp. a^' denote the dynamics with initial condition a 
resp. a' using the same r and U, with the above construction, one has P-a.s. 

Vi>0 a^{t)>a^'{t). (3.2) 

It also yields monotonicity with respect to boundary condition and other nice properties. 

3.2. The height function of the simple exclusion process. We deal in this subsection and 
the next ones, of special initial conditions and boundary conditions, for which the the interface 
between + and — at all times is given by the graph of a function in the coordinate frame (fi, (2) 
(recall the definition above equation (2.17)). This cases are easier to treat for two reasons: firstly 
the rescaled interface motion can be written a functional PDE, which is easier to deal with than a 
flow of curves; secondly, there are bijective correspondence with particle systems that facilitates 
the probabilistic treatment of these interface problems (and also gives extra motivation to study 
them). 

Results concerning interface models are one of the principal building bricks for the proof of 
Theorem 2.2 and 2.4. We describe the correspondence in this subsection and state convergence 
results in the symmetric and the asymmetric cases in Section 3.3 resp 3.4. 
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We say that a set A C (Z*)^ or A C is increasing if and only if 

yxeA,{y>x^yeA), (3.3) 
where the order in (Z*)^ is the usual partial order (2/1,1/2) > {xi, X2) where yi > xi and 2/2 > X2- 
We define the interface between + and — for a configuration a to be the topological boundary 

of 

A-{a):= U (3.4) 

{x: crx = -l} 

which also is equal to the boundary of A~^{a) with an analogous definition. 

If the set of + spins is an increasing set for cr(0) then it remains an increasing set for cr(t) for 
all t > (the only site whose spin can flip have two + neighbors and two — neighbors and the 
reader can check that flipping one of them does not break the property). Also, when the set of 
+ spin is increasing, the interface between + and — is the graph of a function ?/ : R — )■ R in the 
frame (0, fi, £2). The restriction of r/ to Z is integer valued and satisfies \rj(x + 1) — r](x)\ = 1 for 
all 2; G Z and it is linear on every interval of the type [x,x + 1], x £ Hence in a small abuse 
of notation we can also consider 77 as a function on Z, and see the dynamics as a Markov chain 
on the state-space 

0° := {?7 G Z^ I \r]{x + 1) - ??(x)| = 1, Vx G Z} . (3.5) 

In the case where the dynamic is restricted to Ax,, we choose the boundary condition to be + 
on the up and right sides and — on the two opposite side and the interface function rj defined 
on [-2L,2L] (see Figure 2). 
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Figure 2. One-to-one correspondence between the dynamics in a x M 
rectangle with mixed boundary conditions and the corner-flip dynamics on paths. 
This correspondence is of course also valid in the case of a 2L x 2L square. An 
example of possible flip update is displayed 



The dynamic of the interface can be described as follows: sites x G Z (or x G [—2L + 1, 2L — 1]) 
are equipped with independent Poisson clocks with rate one; when a clock rings at x the path rj is 
replaced by rj^ with probability e''/(2 cosh(/i)) and rjx with probability e''/(2 cosh(/i)), where rj^ 
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(and rjx) are respectively the maximum and minimum path in Q that coincides with r] on Z\{x}. 
The paths i]x and r/^ differ if and only if 7] has a local extremum at x i.e. if r]{x + 1) = r]{x — 1)). 

This description is only formal in the full line case since the set of update time is dense, but 
is a rigorous definition of the generator of the chain (for a general proof of existence of a chain 
with such generator see [20, Chapter 1]). 

Finally this interface dynamic can be mapped onto a one dimensional particle system. For 
any a; G Z set 

C(x) = (?7(x) - 7?(x + 1) + l)/2 

and say that a particle lies on the site x when ^(x) = 1. With the alternative view the dynamic 
can be described as follows: each particle jumps to the right with rate e^/ (2 cosh(/i)), and to the 
left with rate e~'^/(2 cosh(/i)), with jumps being canceled if the aimed site is already occupied 
by another particle. In the case of r] defined on a segment, we define the particle in the same 
manner, with the constraint that the extremity of the particle system are "closed", so that 
particle cannot run through them. This system is called the simple exclusion process . 

When h = 0, the particles perform symmetric motions and one speaks about SSEP (symmetric 
dimple exclusion process). When h G (0, oo], the particles are biased towards the right and one 
talks about ASEP (asymmetric). When h = oo particles can only jump to the right and the 
system is called TASEP (totally asymmetric) 

For all these particle systems, the only infinite volume equilibrium measures are the one under 
which ,^(x) are IID Bernoulli variables. 

3.3. Scaling limit for the Asymmetric Simple Exclusion Process. This relation to one 
dimensional particle systems gave additional motivation to study these interface motion and 
result for the scaling limit of the interface have been obtained. Let us cite the first one, due 
to Rost [24] that concerns the case where the initial condition is — in a quadrant and + in the 
three others, in the totally asymmetric case (recall the definition of the function g in equation 
(2.16)). 

Theorem 3.1 ([24, Theorem 1]). Consider the stochastic Ising model on (Z*)^ at temperature 
zero with h = oo and initial condition ctq — in and + elsewhere. Let r/(x, t) denote the 
function whose interface in (0,fi,f2) is the interface between + and — for a{t) (r]{x,0) = \x\). 
Then one has for all e > and T < oo 



lim P 

L— >oo 



max \—r]{Lx, Lt) — g{x,t)\ >e 
te[o,T] 



(3.6) 



The above the result has been developed in order to be able to treat all kind of initial condition. 
Let uq be a 1-Lipshitz function and let (Tq be a sequence of initial configuration for which the 
set of — is an increasing set and the initial interface function t/q satisfies: 



max 



0(1). (3.7) 
Then define u{x, t) as being the unique viscosity solution of 

(3.8) 



'dtu = i(l - (a,u)2), 



u(-,0) =uo, 
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which can be showed to be equal to (recall (2.16)), 

u{x, t) := mf{uo{y) + g{x - y, t)}. (3.9) 

y 

We use (3.9) as the definition of u. It turns out that u is the description of the scaling limit of 
the interface when the initial condition satisfies (3.7). Before stating the theorem, let us explain 
briefly on heuristic grounds why it is so. 

With use the particle system description given at the end of the previous section. With this 
setup, the height variation r]{x, t) —t]{x, 0) is equal to twice the number of particle jumping from 
X — 1 to X in the time interval [0, t]. 

We consider the simplified case where the initial profile is linear with slope s and the particle 
system is in an equilibrium configuration, i.e. {(,x)x£Z are IID Bernoulli variables of parameter 
p := (1 — s)/2. In that case the jump rate of particle from 2: — 1 to x is equal to 

P[r?(x - l,t) = 1 ; r?(x,t) = 0] = p(l - p) = ^(1 - s^). 

Assuming that there is some kind of ergodicity in the system, a law of large number should hold 
and for any fixed x a.s. 

i]it, x) - 7?(0, x) = ]^{l- s^)t{l + 0(1)). (3.10) 

The reason why the argument also works when the original density of particle is non uniform is 
that locally, the system relaxes quickly to equilibrium so that the field ?/(x, t)x^z looks locally 
like IID Bernoulli after a small time. 

The following Theorem was proved in [23] for ASEP in arbitrary dimension. The reader can 
also refer to [25] for a generalization to the /C-exclusion process. 

Theorem 3.2. Consider the stochastic Ising model at temperature zero with h = 00 and ini- 
tial condition ctq as above (and rj^ the corresponding interface). Then, the resettled interface 
converges in law to u in the following sense: For every e, for every positive T and K 

= 0. (3.11) 

These two results remain valid for the ASEP, but time has to be rescaled by a factor 7(/i) = 
coth(/i). 

3.4. Scaling limit of the Symmetric Simple Exclusion Process. The case where h = 
corresponds to the case where the particles perform symmetric jumps. In that case the speed of 
the particle is zero, and one has to rescale time by to get a non-trivial scaling limit. 

It is now a classical result that any dimension, the weak-limit of the density profile of particle 
for the SSEP is given by the heat equation (see e.g. [11, Chapter 4]). In [15] we proved an 
analogous strong limit result for the profile 77, in the case of dynamics restricted to a box: consider 
the zero-temperature dynamics with zero magnetic fields on with + boundary condition on 
the up and right sides and — on the two others. Let 1 1— )• 7]{-, t) defined on [—2L, 2L] denote the 
function whose graph in (0, fi,f2) is the interface between -|- and — for the zero temperature 
dynamics. 

Given a 1-Lipschitz function vq : [—2,2] 1— )■ M with t;o(±2) = 0, assume that one start the 
dynamics in A/, with a sequence of initial condition (Jq for which the interface rj satisfies 



lim P max [—77 {Lx,Lt) 

L-s-oo \x\<K L 
. t<T 



u{x, t)\ > e 
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lim sup 

L^'^ x&[-2L,2L] 



Let V be the solution of 




-ri{x,0) - vo{Lx] 



vo, VxG[-2,2], 
0, Vt > 0. 



(3.12) 



(3.13) 



Theorem 3.3 ([15, Theorem 3.2]). Consider the dynamics on with the above mentionned 
initial condition, then for all T >0 and e > 0, w.h.p. 



sup 

t<=[0,T],xe[-2,2] 



—ri{Lx, L^t) — v{x, t) 



< e 



(3.14) 



where v is the solution of (3.13). 



Let us give shortly a heuristic explanation for Theorem 3.3. The first thing is to show that 
the expected value of (??x)x-e[-L,L] satisfies approximately the heat-equation. 

First notice that the expected drift of r]x at time zero depends only on the value of r]x and 
r]x±i. If rjx is a local maximum it will jump down by two with rate 1/2 whereas if it is a local 
minimum it will jump up by two with the same rate. Else r]x has drift zero so that 

1 i]x±i{0) = r]x-i{0), 

dtnVx{t)]\t=o = 1-1 if Vx±i{0) = r/,+i(0), (3.15) 
else. 

v 

The reader can check that the r.h.s. of the above equation is equal to 

^{Vx+i{0) + Vx-i{0) - 2r/,(0)) =: ^A^rj^, 

where denotes the discrete Laplacian. Hence, using the Markov property, one obtains that 
mr]x])x(^[-L,L] satisfies 



dtnr]x{t)] = ^AdE[7?,.(t)]. 



(3.16) 



Theorem 3.3 is obtained then by showing that: 

• The solution of the discrete heat-equation converges to the solution of the continuous 
one in the scaling limit (this is classic). 

• That r]x{t) concentrates around its mean for large values of L. This is more delicate. We 
proved it by proving concentration for all the Fourier coefficient in a base of eigenfunc- 
tions of Afj. 

Projections and trigonometry allows a heuristic derivation of (2.7) from Theorem 3.3 and 
equation (3.13). Indeed, it is quite reasonable to think that for the original dynamics with 
shrinking — domain, the local drift of the interface is the same than for this dynamics with 
modified boundary condition. However, there is a crucial ingredient missing to try to perform a 
proof. 

Indeed Theorem 3.3 does not say anything about the drift of the interface around the "poles" 
of Aiit), i-e. around the points for which one of the coordinates is extremal. This problem was 
treated using another correspondence with particle system, first in [27] where a complete sketch 
of proof was given in a special setup (periodic boundary condition). This study was pushed 
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further in [15] where all the technical details were handled to get a result that could be used to 
prove Theorem 2.2. We present this correspondence in the next Section. 

3.5. Dynamics near the "poles": Zero Range Process and Scaling Limit. Near the 
"poles", the dynamics cannot be easily reduced to an interface dynamics: indeed (see Figure 
3) there are some possibilities for the set of — to break into several connected components. 
However one can introduce an auxiliary dynamics that cancels the transitions that makes AL{i) 
disconnected. 
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Figure 3. An example of spin update that splits the interface into two discon- 
nected components. The interface dynamics presented in this section does not 
allow this kind of move. 



More precisely, we consider with boundary condition — in the upper-half plane and — 
one the lower-half plane, staring with an initial condition such that the interface between -|- and 
— is the graph of a function [—L,L] — ?• M (plus some vertical lines), and run a modified Ising 
dynamics that discards update if the interface after the update is not a single connected curve. 
With this dynamics, the interface remains the graph of a function for all time (if one neglects 
the vertical lines). We call r](x,t) the corresponding function (as there is no confusion possible 
with r] from the other section): by convention we choose it to be defined on [—L, L] D Z* as it is 
piecewise constant. 

In this case also, we can describe the evolution of the gradient as a particle system. For 
X G {—L, . . . ,L} we set 

Cx{t) ■■= Vx+i - Vx-l 

to be the discrete gradient of t]. We say that each site in {—L, . . . , L} carries |'^x(i)| particles. 
These particles are said to be of type ^ if ^ is positive and of type S if ^ is negative. 

Under the modified Ising dynamics depicted above, the rules for the motion of the particles 
is the following: 

• If there are k particles on a site, they jump left or right with rate 1/(2A:). 

• If a particle of type A meets a particle of type B they annihilate. 

This kind of particle system where the jump rate depends on the number of particles on one 
site is called zero-range process (see Figure 4 for a scheme of the correspondence of interface 
dynamics with the particle system), and has been extensively studied in the literature (see e.g. 
[1] where invariant measure of the process are studied). 

This correspondence with particle system was underlined in [27], and a partial proof of the 
scaling limit of the interface motion was given there: the scaling limit of ^ should be the solution 
of the equation 
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Figure 4. Correspondence between interface dynamics and zero-range process. 
Arrows represent possible motions for the interface and their representation in 
terms of particle moves. When an A particle jumps on a B particle (green arrow), 
they annihilate. 



dtw = -^^^. (3.17) 

Although we were not able to complete this proof fully (in particular we miss a statement 
concerning existence and regularity of the solution in (3.17)), we proved a partial statement 
that was sufficient for the purpose of the proof of Theorem 2.2. Our statement is that, on the 
macroscopic scale, r]{x, t) stays close to a deterministic discrete evolution that can be though as 
a discretization of (3.17) [15, Theorem 3.4]. We record informally here the result for quotes in 
the rest of the paper 

Theorem 3.4 ([15, Theorem 3.4 and Corollary 3.5]). In some weak sense equation (3.17) 
describes the evolution of the rescaled interface on the diffusive time-scale. 

4. Zero magnetic fields: a detailed sketgh of proof for Theorem 2.2 

We expose in this section the ideas behind the proof of Theorem 2.2. The first important 
step is to reduce the proof to an infinitesimal statement. Using continuity properties of the 
conjectured scaling limit, we can show that in order to control the evolution of the + domain 
for all time, it is sufficient to control the motion with a first order precision during a small time 
e. 

Once this is done, it suffices to iterate the statement many as much as need (order e~^) to 
control the evolution for arbitrary positive time. We do not develop this point further. 

Let us state directly the two infinitesimal statement we want to prove: the first one concerns 
continuity of the interface motion (which has to be used to get continuity of the motion). 

Proposition 4.1. Let T> be convex and with a Lipschitz curvature function. For every a > 0, 
w.h.p. (recall definition (2.8) j 

AiiL'^t) C LV'-''^ for every t>0. (4.1) 
Moreover, for every a > there exists ei{a, A;max) > such that w.h.p 

AULh) D for every te [0,ei]. (4.2) 



the second one is the control at first order of the evolution. 
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Proposition 4.2. For all 6 > there exists eo(5, ^mim ^max) > such that for all < e < Eq, 
w.h.p., 

AiiL^e) C LV{e{l - 5)), (4.3) 

and 

AiiLh) D LV{e{l + 6)) (4.4) 

where T){t) is the solution of (2.7). 

The main work that remains is to prove the two inclusion bounds (4.3) and (4.4), as Propo- 
sition 4.1 is more of a teclmical detail. We detail now the sketch for the proof of the upper 
inclusion (4.3). The other bound is proved similarly but for technical reason it is a bit harder 
to expose its proof. 

Sketch of the proof of (4.3). What we have to show is that after a time e all the spins in T> \ 
'D(e(l — 6)) (on the rescaled picture) that were initially — have turned + w.h.p. after a time e. 

To do so we divide T) \ 2?(e(l — 6)) in eight regions and in each of these region, we will try 
to compare our dynamics with some interface dynamics. We consider four regions around the 
poles named {Ai)f^^ and four others named {Bi)f^^ (see Figure 5) and one wants to show that 
each one of them is filled with + after a time e. In what follows we will chose the Ai to be very 
small zones around the poles whereas the Bi will cover a proportion of the boundary of P close 
to 1. 

Due to rotational symmetries of the problem, the inclusion (4.3) reduces to proving that 
w.h.p. 

Vx G Bi,ax{£L^) = +1, 

o (4.5) 
yx(^Ai,a^{eL^) = +l. 

The idea behind this division is that initially, 7 restricted to the region Ai (resp Bi) is the 
graph of a function in the coordinate frame (61,62) for Ai and in the frame, resp. (fi,f2). We 
want to modify our dynamics a little bit to get in the context of Theorem 3.3 or Theorem 3.4. 

The idea to prove each line of (4.5) is to replace the dynamic a by another one that creates 
more — spins but that we can handle better. 

Let us start with the case of zone Bi. We look at the dynamics restricted to the rectangle of 
Figure 6, and we decide to cancel all the updates that concerns spins of the down or right side 
of the rectangle, and remark that the spins on the two opposite sides remain + for all time. The 
obtained dynamics a^^^ has more — spins that the original one (it can be seen from the graphical 
construction of Section 3.1) and falls in the setup described in Section 3.4. More precisely we 
have a rectangle with mixed boundary condition instead of a square, but Theorem 3.3 also apply 
in that case). 

If Tj^^^ denotes the interface function corresponding to a^^^ (in the setup of Section 3.4) and 
q{-,t) denote the function whose graph in the coordinate frame (fi,f2) is the boundary of I'(t) 
restricted to the zone Bi (we chose t small enough so that q remains defined on some "large 
interval" /). The first line of (4.5) is proved if one can show that 

\r]^^\Lx,L^e) < q{x, £{1 + 5)),yx G /, (4.6) 

where / is an interval that is strictly included in the domain of definition Iq of q{-,t) which 
sufficiently large to have control on the whole region Bi. 
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Figure 5. Decomposition of P \ T>{e{l — S)) in eight zone, four small ones 
around the poles Pi that we call Ai and four larger one away from the poles, 
called Bi 



Thus using Theorem 3.3, it is sufficient to prove that 

u{x,e) < q{x,e{l + 5)),yx e I, (4.7) 

where u is the solution of equation (3.13) with initial condition uq = q{-,0) corresponding to the 
initial position of the interface. Continuity properties of the heat equation allows us to say that 
if the interval I is fixed and uq is smooth enough, one has uniformly in x G /, 

u{x,e) = no(x) + ^d'^uo{x)e + O(e^). (4.8) 

This does not hold uniformly in the full interval Iq because of the Dirichlet boundary condition 
that makes the drift equal to zero at the extremities of the interval. The anisotropic motion by 
curvature is sufficiently regular to obtain something similar for q{x,e): uniformly on x € / 

q{x,t)=t^dlq{x,0)+O{t^). (4.9) 

At an informal level this is just a Taylor formula combined with (2.7) and some trigonometry 
(for the projections). 

As curvature remains positive everywhere for all time and q{x,0) = u{x,0), it is straight- 
forward to deduce (4.7) from (4.8) and (4.9) for small enough e. 

The treatment of zone a around the pole (say Ai) is more delicate but follows the same line. 
The difference is that we have to perform a longer chain of modification of the dynamics. Each 
modification makes the dynamics have more — spins, or result in a dynamics that coincide with 
the previous one with large probability. 
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Spins that remains + for all time 




Figure 6. In order to control the dynamics on a zone of type A, we look at 
the dynamic restricted to a "quadrant" of the original shape. In order to cut 
dependence from what happens outside the rectangle, we block the update on a 
set of spins that have to stay — at all time (dark red) , and we can do a monotone 
coupling of this dynamics with the original one. On the two opposite sides of the 
rectangle (dark blue), even if updates are not rejected, note that the spins have 
to stay + for all time due to the majority rule. 



We look at a dynamic restricted to a rectangle around the pole that includes Ai, e.g. take 
the rectangle twice as large as Ai (see Figure 7) . We want to modify the dynamics by fixing the 
boundary condition on the rectangle (or to freeze updates) . The problem here is that to get to 
the setup of interface dynamics of Section 3.5 we need to freeze some sites with + spins as well 
as sites with — spins so that the modified dynamics does not compare well with the original one. 

A way to get around this problem is to prove first that sites at a macroscopic distance (i.e. 
positive distance on the rescaled picture) from T>{0) do not change sign ever with large probability 
(this is for instance the upper inclusion bound in Proposition 4.1). 

Knowing this we can freeze spins to + in the upper-part of our rectangle provided it does 
not touch I?(0): the obtained modified the dynamics will coincide with the original with large 
probability. 

Then if we freeze the spins in the lower-part of the rectangle to — we get a dynamics that 
compares well with the previous one (it has more — spins, see upper part of Figure 7 upper 
part). 
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Once the whole boundary has been frozen, we add — spms to our initial condition so that on 
the macroscopic scale, our initial condition is a smooth interface: we cancel the irregularity at 
the boundary (see figure 7 lower part). 

The dynamics restricted to the rectangle is not exactly an interface dynamics so we have 
to modify it once again. We cancel all the moves that break the interface between + and — 
into several connected components. With our modified initial condition (an interface that has a 
unique local maximum) a disconnection of the interface can only occur by adding a + spins so 
the modified dynamics has one again more — spins than the original one. 

Let (T^^^ denote the last mentioned modified dynamics and i]^'^\x,t) denote the interface 
function that corresponds to it. Let r denote the function whose graph in (ei, 62) is the boundary 
of T)[t) in the zone Ai. 

We are left with proving that 

yr/(2)(Lx,L2e) < r(x, e(l + 5)), Vx G J, (4.10) 
ij 

for a small interval J around the pole. 

Using Theorem 3.4, one gets that for small a > 0, for any positive t w.h.p. 

y7?(2) (Lx, L^t) < w{x, (1 + a)t), (4.11) 

where w is the solution of dtw = {l/2)d'^w with Dirichlet boundary condition, and initial 
condition given by the the interface. The number a can be taken arbitrarily small by reducing 
the size of the zone around the pole that is considered. 

With this information, proving (4.10) reduces to show that 

w{x, (1 + a)E) < r{x, e{l + 8)) Vx G J, (4.12) 

where / is an interval corresponding to the zone Ai. Equation (4.12) is proved in the same 
manner that we proved (4.7), we choose a = (5/2, and perform Taylor expansion at first order 
on both side. 



w{x, (1 + 5/2)e) = w{x, 0) + e(l + 5 /2)/2{d^fw{x, 0) + 0{e^), 

Again the second line is formally (2.7) and some trigonometry. Then using the fact that w{-, 0) = 
r(-,0) on an interval around the pole and that |(93;^(a^; 0)| is uniformly small on J if J is chosen 
small enough, we can deduce (4.12) for (4.13) for e small enough, with a small interval J. 

Proving the other inclusion relies on the same ingredients (4.4): first separate different zones 
that one wants to control, and then, for each zone reduce to an interface dynamics via a chain 
of modification. After doing this we use the scaling results Theorem 3.3 and 3.4 to conclude. 
Then one derives the result for the original dynamics using monotonicity. The chain of dynamic 
modification is a bit longer and tedious in this case that for the upper-bound (4.3). 

□ 
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Spins that remains + for all time with large probability 
i_ 

Spins for which updates are blocked (stay — /| ^ 




adding — spins to get $. smooth boundary condition on the rescaled picture 




Figure 7. In order to control the dynamics on a zone of type B, we look at the 
dynamic restricted to a rectangle around the pole. In order to cut dependence 
from what happens outside the rectangle, we block the update on a set of spins 
that have to stay — at all time (dark red). On the rest of the boundary we note 
that the spins have to stay + for all time with large probability so that we can 
couple the dynamics that reject updates on those site so that it coincides with 
the original dynamics with large probability. 



5. Positive magnetic fields: proof of Theorem 2.4 

The proof of the scaling limit of the — droplet in this case uses several techniques in common 
with the zero-magnetic field case, but several features of the model makes it much simpler in a 
sense: 

• The scaling limit, although being non-smooth, is a much simpler object than anisotropic 
motion by curvature. In particular there is no instantaneous travel of information on the 
rescaled picture: if one modifies one side of the droplet, it will take a positive time to 
have an effect on regions that are at a positive distance of where the modification took 
place. 

• + spins remain + for all times, making the set of — a decreasing set in time. 

Recall that the aim of this Section is to prove Theorem 2.4 in the special case where cj(0) is 
equal to - in [-L,L]2 n (Z*)^ and + outside of it. Recah that in that case T>{t) = nf^iA(t) 
where Pi(t) is defined in equation (2.17). 

The upper inclusion j^Aiit) C V^^^t) is an easy consequence of Theorem 3.1: the idea is to 
couple the dynamics with four corner growth dynamics using the same clock process for updates. 
What can be rather surprising is that this simple method gives a sharp bound. 

Indeed for the model based on simple exclusion, one can get an upper bound in the same 
manner by comparing with four corner growth dynamics, replacing g[x, t) by the solution of the 
heat equation with initial condition |x| (and time scaling L by L?), but this would not be sharp. 
Indeed the shape obtained as an upper-bound with such a technique presents some angles on 
points with extremal coordinate (poles), whereas the true scaling limit is smooth (at least C^). 
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The reason is that the drift of the interface at the pole is positive when the pole is convex, and 
one needs to apply Theorem 3.4 to take the drift into account. 

Here on the contrary, nothing happens around the poles due to the singularities of the function 
b{d) (recall (2.12)) that is equal to zero for = i7r/2, i = 1, . . . ,4. The main thing to prove to get 
the lower bound j^Atii) ^ 'D^~^\t) is to show that the interaction between the four quadrants 
of our shape is quite limited. 

It is quite difficult to control directly what happens around the pole when the interface is not 
completely flat, but one finds a way to bypass this problem. Similarly to what is done in [14] 
in the case of dynamics for polymer with an attractive substrate, one adds, in a quite artificial 
manner, a fiat part of interfacer around the pole, and modify a bit the statement that has to be 
proved (see Proposition 5.4 below). Adding this fiat part makes the evolution of the four corners 
almost independent of one another for some positive time, and thus allows to use Theorem 3.2 
to control the evolution of each corner and thus of the total shape. 

5.1. The upper-bound. The upper-inclusion of Theorem 2.4 follows quite easily from Theorem 
3.1. We prove it in this section. 

Proposition 5.1. For any 6 > 0, for the dynamic with h = oo starting from the initial condition 
— in [— L,L]^n (Z*)^ and + outside of it, one has w.h.p. 

{jALiLt))t>o C P(^) (t), Vt > 0, (5.1) 
where V{t) is defined by (2.17)-(2.18). 

Proof. Using the graphical construction of Section 3.1, we can couple the dynamics a{t) with ini- 
tial condition a^: — in [—L, L]^ (we drop intersection with (Z*)^ in the notation for conciseness) 
and + elsewhere with other dynamics using the same clock process: 

• The dynamics ai with initial condition — in [— L,oo)^ and + elsewhere. 

• The dynamics 1T2 with initial condition — in [— L,oo) x [— oo,L) and + elsewhere. 

• The dynamics 1T3 with initial condition — in [— 00,-L)^ and + elsewhere. 

• The dynamics 1T4 with initial condition — in [— oo,L) x [—L,oo) and + elsewhere. 

We define A\{t), A?]^{t), ^|(t), A\{t) analogously to Aiit) of (2.3) for these four dynamics. 
According to monotonicity properties of the dynamics in the initial condition one has AL{t) C 
A\{t) for all i e [1, 4] so that: 

AL{t) c Al{t) n Al{t) n Al{t) n Alit). (5.2) 

The dynamics is the same as the one considered in Theorem 3.1 with a space shift of the 
initial condition, and thus by Theorem 3.1, 

Al{t) c {xfi + yh\y> g{x,t) -2-5}=-. Ai{t). (5.3) 

Defining Af{t), i = 2,3,4 as rotations of Af{t) by angle {i— l)7r/2, one gets analogous inclusion 
for A''^{t) by symetry, and using (5.2) we get that w.h.p. 



AUt) c Af{t) n Ai{t) n A^^{t) n 4(t) c v^^\t). 



(5.4) 

□ 
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5.2. The lower bound. The aim of this section is to prove the other inclusion of Theorem 2.4 

Proposition 5.2. For any 6 > 0, for the dynamic with h = oo starting from the initial condition 
— in [— L,L]^n (Z*)^ and + outside of it, one has w.h.p. 

{^AL{Lt))t>o D V^'-^^t), Vt > 0, (5.5) 
where V{t) is defined by (2.17)-(2.18). 

Let us explain how the proof goes. Consider the sets A\^{t) defined in the previous Section. 
The idea in our proof of the lower bound it to show that the inclusion used in (5.2) is almost 
an equality. We will decompose the proof in two steps: 

• First, we prove that when the rescaled-time t is smaller than one the inclusion (5.2) is 
indeed an equality with large probability. 

• Second, when the rescaled time t is larger than one, we use a special strategy involving 
adding portions of straight-line on the interface around the pole in order to reduce oneself 
to a case where one can treat the four corner dynamics independently. 

The control of the dynamics for times t < 1 is summarized in the following lemma 

Lemma 5.3. For any 5 < one has with high probability, for all t < 1 — 6 one has 

Vl^'^ C jAULt). (5.6) 

Proof. Define ICl to be the sets of sites in [— n (Z*)^ that are neighboring either the 
horizontal or vertical axis, 

/Cl := {x G (Z*)2n I min(|xi|,|x2|) = 1/2}. (5.7) 

Set 

r := inf{t >0\3x£ lCL,CF^{Lt) = +}. (5.8) 
The reader can check that for t < t, the evolutions of the four corner do not interact with one 
another and that 

AUt) = Alit) n Alit) n Ai{t) n AUt). (5.9) 
Moreover, as a consequence, one has 

r := inf Tj := inf{i > I 3x e /CL,3i e [1,4], o-j,(Lt) = +}. (5.10) 

i=l..A 

What remains to check is that ti is roughly equal to L. If one uses the particle system 
description from Section 3.2, the dynamic cJi corresponds to TASEP with step initial condition: 
i.e. with the negative half-line full of particles and the positive half-line empty. 

With this description, ti is the first time that either the rightmost particle, which starts from 
0, reaches L, or the leftmost empty space (or antiparticle) hits —L. As the jump rate of the 
particles (and thus of antiparticles) is equal to one, the central limit Theorem for simple random 
walk gives us that these time are equal to L + 0(L^/^) where the second term includes the 
random normal correction. Thus one can infer that w.h.p: 

n>L-L'/\ 

so that the same is true for r. 

From Theorem 3.1 one gets that with high probability for all t < 1, 



Al{Lt)D{xii+yi2 \ y>g{x,t)-2 + 6} = A^^{t). 



(5.11) 
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SO that combined with (5.9) one gets that w.h.p. for all f < 1 — 5, 

AULt) D A-^\t) n A^\t) n A^\t) n A-^\t) d v^^^\t). (5.12) 

□ 

We address now the issue of time larger than one. At this time, the different corners start to 
interact so that one has to find a trick to regain independence. For t S [1, 4], set d{t) = 2\/t — t. 
The reader can check that for all t > 1, T>{t) is inscribed in the square [—d{t), d{t)Y ■, i.e. 

max X • ei = — min x • ei = dit), (5.13) 

x6l?(t) xgl'(t) 

the same being valid for the second coordinate. 

We will intersect T){t) with a smaller square, which gets us a shape with vertical and horizontal 
edges around the poles. More precisely, given 5, for t < 4(1 — 6) (we have in that case d{t) > 6 
if 6 is small enough) we define the function 

.^y^[9ix,t)-2ii\x\<dit)-5, 

\\x\-{d{t)-6)+g{d{t)-6,t)-2ii\x\>d{t)-6. ^' ' 

It is continuous in x, convex in x and has slope equal to ±1 outside of [—d{t) + 5, d{t) — 6]. Define 
Pi(t) to be the epigraph of g{x, t) in (0, fi, (2) and Pj(t) its rotation by an angle {i — l)7r/2, and 
set 

4 

P(f) = f|A(t). (5.15) 

i=l 

The reader can check that 

V{t)=V{t)n[-r{t),r{t)]\ (5.16) 
where r(t) G id{t) — 6/2, d{t)) is the unique solution of the equations 

g{x, t) = —X. (5-17) 
To check that r(t) > d{t) — 6/2 it is sufficient to observe that 

g{d{t) - 6/2, t) = 5/2 + g{d{t) - 6,t) = 5/2 - d{t) + "MO'^ + ^ ^^^^ _ ^^2. (5.I8) 

Hence P(t) has flat parts around the pole. Thus, if one starts a dynamics from a shape 
approximating T>{t) the four corner will not interact instantaneously like in the proof of Lemma 
5.3. We will use this fact to prove the following result 

Lemma 5.4. For any 5 < and e small enough one has for all k > such that ek < 3 — 45 
w.h.p, 

yAlHI -5){1 + ek)L) D V{1 + ek). (5.19) 
As a consequence, for all t>\ — 6 one has 

\AL{Lt)z>V{t)^-^^\ (5.20) 
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d{t) -r{t) 



r{t) 



d{t) 



Figure 8. The shape T>{t) (darkened on the figure) is obtained by intersection 
four rotated version of the epigraph of the function x i— t- g{x, t) — 2 (recah (2.16)). 
For t G [1,4] (as above) this shape is inscribed in [—d{t),d(t)]'^ where d{t) = 
2^/t-t. We define V{t) to be the intersection oiV{t) with [-r{t),r{t)f defined 
in (5.17). 



Proof. The main part of the job is proving (5.19), as (5.20) is deduced from it by monotonicity 
in t of ^L(t). We proceed by induction on k. For /c = the result is just a consequence of 
Lemma 5.3 for t = \ — 6. 

Now we suppose that the result is true for k and and prove it for k + 1. If (5.19) holds for 
k w.h.p, then, using the graphical construction, one can couple ((t((1 — 5)(1 + ek) + t))t>o with 
a^{t) a dynamic starting with initial condition — in LV{1 +£k) and + elsewhere, in such a way 
that 



AL{{l-5){l + ek) + t)Z) AUt). Vt>0, 



(5.21) 



whenever 



a{{l + ek)) D LV{1 + ek), 

where A\ is defined as in (2.3) with a replaced by . 

As (5.22) holds w.h.p. we can prove (5.19) for A: + 1 if one proves that w.h.p. 



(5.22) 



Al{{l-5)e)Zi LV{l + e{k + l)) 



(5.23) 



Now, we couple (J^{t) with four interface dynamics a^'^'^\t) starting with initial condition with 
initial condition — in LT>i{l + ek) and + elsewhere using the same clock process. 

As we remarked in the proof of Lemma 5.3, one has 



4 




(5.24) 



i=l 
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and this inclusion is an equality up to the first time a site near the axes changes spin. More 
precisely set 

Kl := {x G n [-Lr{t), Lr{t)f \ min(|2;i|, jxal) = 1/2} = /Cl n LV{1 + ek), (5.25) 

and 

f := mf{t >0\3x£ /Cl, a^XO = +}• (5-26) 
One has from the definition of our dynamics that 

Ai{t) = f]Ai''\t), yt<f. (5.27) 
1=1 

Our first task is to show that w.h.p. f > e in order to be allowed to use (5.27). We check that 

fi := {t > I 3x G iCL,(J^^'^\Lt) = -}>£, (5.28) 
which by symmetry is sufficient. Indeed 

f = min = Tj 

iG[l,4] 

where the fi are defined similarly to fi. 

We use again the correspondence with particle system of Section 3.2 to do that: the dynamics 
corresponds to a particle system with an initial condition where the right-most particle is 
located at L{d{l + ek) — 5) and the left-most empty-space at —L{d{l + ek) — 5). 

The time fi is the first time where either the rightmost particle hits Lr{l + ek) or the leftmost 
antiparticle hits —Lr{l + ek). Using central limit Theorem for the sums of exponential variables 
one has w.h.p, 

fi>L{{r-d){l + ek) + 5) + L^/^. (5.29) 
As r(t) > d{t) — 6/2, one has fi > S/AL > eL w.h.p. provided that e < (5/4. 

Now, having shown that equation (5.27) holds for t = e w.h.p., to get (5.23) it is sufficient to 
prove that w.h.p. (recall the definition of T>{t) in (5.15)) 

A^I'''\{1 - 5)eL) D Pi(l + (A; + l)e). (5.30) 

As the two sets are epigraphs of function in (0,fi,f2) it is sufficient to prove an inequality 
between function. Let rj^{x,t) denote the interface function corresponding to cr^^'^). Theorem 
3.2 give us the scaling limit for the evolution of r/^(x, t), which is given by u{x, t) the solution of 
(3.8) with initial condition 

uo{x) := g{x,l + ek). (5.31) 
Thus (5.30) is proved if we can show that 

u{x,{l-6)e) <g{x,l+e{k + l)), G M. (5.32) 
The reader can readily check, using (3.9) that 

u{x, (1 - 6)e) = g{x, 1 + e{k + 1) - 6e) - 2 < g{x, 1 + e{k + 1)), 

Vx G [-d{l + ek) + 6, d{l + ek) - 6]. (5.33) 
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What remains to show is that the inequahty is also vaUd outside of the interval [— d(l + ek) + 
6,d{l + ek) — 6]. Set x > d{l + ek) — 5 (by symmetry of the functions it is sufficient to check 
this case), as u{-,t) is 1-Lipshitz for all t 

u{x, (1 - 6)e) < u{d{l + ek) - 5, (1 - 6)e) + (x - d(l + ek) + 6) 

< g{d{l + ek)-6,l+ e{k + 1)) + {x - d{l + ek) + 5)= g{x, 1 + e{k + 1)), (5.34) 

where the last inequality comes from the definition of g. 

Outside of this interval the inequality (5.32) is also valid: indeed in that region u(x, (1 — 6)e) 
is 1-Lipshitz in x and g{x, l + e{k + l)) has slope +1 and —1 on the right resp. left of this interval 
(5.33). This concludes the proof of (5.30) and thus of (5.19). 

We can now detail the proof of (5.20). For t > 4 - we note that that T>^~^^\t) = so that 
the inclusion is trivial. 

For t G (1 — (5, 4 — 45], note that, with high probability, (5.19) holds for ah k with efc < 3 — 46. 
Let kt be the smallest integer such that (1 — 5){1 + ek) > t. As AL^t) deacreases in time, from 
(5.19), for alH e (1 - 5, 4 - 45] one has w.h.p. 

^AULt) D ^Al{{1 -S){1 + ekt)) D V{1 + eh). (5.35) 

Then we note that P(t) D V'^'^^t) and that 

l + ekt< + e, (5.36) 

1 — 

so that (5.35) implies 

jAULt) D (^-^ + e^D V^~^^\ (5.37) 

if e and 5 are chosen small enough. 

□ 

Proposition 5.2 is just the concatenation of Lemma 5.3 with Lemma 5.4. 

Acknowledgments: The author would like to thank the organizers of the 2012 PASI con- 
ference, where he had stimulating discussion with other participants, Milton Jara, for valuable 
bibliographic help concerning TASEP, and Frangois Simenhaus and Fabio Toninelli for numer- 
ous enlightening discussion on the subject. This work was partially written during the authors 
stay in Instituto de Matematica Pura e Applicada, he acknowledges the kind hospitality and 
the support of CNPq. 

References 

[1] E. D. Andjel, Invariant measure for the zero-range process, Ann. Probab. 10 (1982), 525-547. 

[2] P. Caputo, F. Martinelli, F. Simenhaus and F. L. Toninelli, "Zero" temperature stochastic 3D Ising model 

and Dimer covering fluctuation: a first step towards mean curvature motion, Comm. Pure Appl. Math. 64 

(2011), 778-831. 

[3] R. Cerf, S. Louhichi, The initial drift of a 2D droplet at zero temperature, Probab. Theory Relat. Fields 137 
(2007), 379-428. 

[4] L. Chayes, R. H. Schonmann and G. Swindle, Lifshitz' Law for the Volume of a 2-Dimensional Droplet at 

Zero Temperature, J. Stat. Phys. 79 (1995) 821-831 
[5] A. De Masi, E. Presutti, E. Scacciatelli, The weakly assymetric simple exclusion process, Ann. Inst. H. 

Poincare Proba. Stat 25 (1989), 1-35. 
[6] I. Dolcetta and P. Lions, Viscosity Solutions and Applications, Springer, 1995. 



28 



H. LACOIN 



[7] H. G. Eggleston, Convexity, Cambridge University Press, New York 1958 

[8] M. E. Gage and Y. Li, Evolving plane curvature in relative geometries II, Duke Math. J. 75 (1994) 79-98. 
[9] M. A. Grayson, The heat equation shrinks embedded plane curves to round points, J. Differential Geometry 
26 (1986), 285-314. 

[10] J. Gartner, Convergence towards Burgers equation and propagation of chaos for weakly asymmetric simple 

exclusion process, Stoch. Proc. Appl. 27 233-260. 
[11] C. Kipnis and C. Landim, Scaling limits of interacting particle systems. Springer- Verlag, Berlin, 1999. 
[12] C. Kipnis S. Olla, and S.R.S. Varashan Hydrodynamic limits and large deviation for simple exclusion process. 

Comm. Pure Appl. Math. 42 (1989) 115-137. 
[13] H. Lacoin, Approximate Lifshitz law for the zero-temperature stochastic Ising model in any dimension, 

(preprint) arXiv:1102.3466 [math-ph]. 
[14] H. Lacoin, The scaling limit of polymer dynamics in the pinned phase, (preprint) arXiv: 1204- 1253 [math-ph] 
[15] H. Lacoin, F. Simenhaus and F.L. Toninelli, Zero-temperature stochastic Ising model in two Dimension and 

anisotropic curve-shortening flow, (preprint) arXiv:1112.3160 [math-ph] 
[16] D. Levin, Y. Peres and E. Wilmer, Markov Chains and Mixing Times, American Mathematical Society, 

Providence, RI, (2009). 

[17] L M. Lifshitz, Kinetics of ordering during second order phase transitions, Soviet Physics JETP 15 (1962), 
939-942. 

[18] E. Lubetzky and A. Sly, Cut-off for the Ising model on the lattice, to appear on Invent. Math. 
[19] E. Lubetzky and A. Sly, Critical Ising on the square lattice mixes in polynomial time, to appear in Comm. 
Math. Phys. 

[20] T. M. Liggett, Interacting Particle Systems, Springer- Verlag, New York, 1985 

[21] E. Lubetzky, F. Martinelli, A. Sly and F. L. Toninelli, Quasi-polynomial mixing of the 2D stochastic Ising 

model with "plus" boundary up to criticality, to appear on J. Eur. Math. Soc. 
[22] L. Onsager, Crystal statistics. I, A two-dimensional model with an order- disorder transition, Phys. Rev. 65 

(1944) 117-149. 

[23] F. Rezakhanlou, Hydrodynamic limit for attractive particle systems on 11^ . Comm. Math. Phys. 140 417-448. 
[24] H. Rost, Nonequilibrium behaviour of a many particle process: density profile and local equilibria, Z. Wahrsch. 

Verw. Gebiete 58 (1981), 41-53. 
[25] T. Seppalainen Existence of Hydrodynamics for the Totally Asymetric Simple K-Exclusion Process, Ann. 

Probab. 27 (1999) 361-415. 

[26] R. H. Schonmann, Slow droplet- driven relaxation of stochastic Ising models in the vicinity of the phase 

coexistence region. Comm. Math. Phys. 161 (1994) 1-49. 
[27] H. Spohn, Interface motion in models with stochastic dynamics, J. Stat. Phys. 71 (1993), 1081-1132. 
[28] N. Sugimine, A lower bound on the spectral gap of the 3-dimensional stochastic Ising models, J. Math. Kyoto 

Univ. 42 (2002), 751-788. 

H. Lacoin, CEREMADE - UMR CNRS 7534 - Universite Paris Dauphine, Place du Marechal de 
Lattre de Tassigny, 75775 CEDEX-16 Paris, France. 
E-MAIL: lacoinOceremade . dauphine . f r 



